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' Abstract 



We present a complete treatment in the strong field limit of gravitational retro-lensing by a 
g3 ' static spherically symmetric compact object having a photon sphere. The results are compared 

, with those corresponding to ordinary lensing in similar strong field situations. As examples of 

application of the formalism, a supermassive black hole at the galactic center and a stellar mass 
black hole in the galactic halo are studied as retro-lenses, in both cases using Schwarzschild and 
Reissner-Nordstrom geometries. 
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1 Introduction 

m 

■ Gravitational lensing is an important tool in astrophysics. Observable phenomena of lensing by stars 
O . and galaxies can be explained in the weak field approximation, that is to say, keeping only the first 

order term in the expansion of the deflection angle 1 . But when the lens is a black hole, a full general 

bfx' relativistic treatment is in order. 
OJ). 

Recently, several articles studying strong field lensing appeared in the literature. Virbhadra & Ellis 
^ ' [2] numerically analyzed the situation where the lens is a Schwarzschild black hole placed in the center 

■ of the Galaxy. They obtained the lens equation with an asymptotically flat metric and found that, 
besides the primary and secondary images, two infinite sets of faint relativistic images are formed by 
photons that make complete turns (in both directions of rotation) around the black hole before reach- 
ing the observer. Fritelli, Kling &l Newman [3j found an exact lens equation without any reference to a 
background metric and compared their results with those of Virbhadra & Ellis. Bozza et al [3] devel- 
oped an approximate analytical method (generally called, as well, the strong field limit) by which they 
obtained the positions and magnifications of the relativistic images for the Schwarzschild black hole. 
Subsequently, Eiroa, Romero & Torres ^ applied the strong field limit to the Reissner-Nordstrom 
black hole and discussed the possibility of detection of the relativistic images in the next few years. 
Bozza jH] extended the strong field limit to any static spherically symmetric lens and analyzed the 
case of a charged Brans-Dicke black hole. Bhadra [7j studied the Gibbons-Maeda-Grafinkle-Horowitz- 
Strominger charged black hole of string theory in the strong field limit. Virbhadra & Ellis numerically 
investigated the lensing by naked singularities [S]. Bozza & Mancini [2! used the strong field limit to 
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study the time delay between different relativistic images, showing that different types of black holes 
are characterized by different time delays. Spinning black holes are more difficult to tackle. Bozza 
[lOj analyzed the case of quasi equatorial lensing by rotating black holes and Vazquez & Esteban 
studied the Kerr black hole where the observer and source have arbitrary inclinations. Many of these 
recent advances (if not all) deal with theoretical implications of General Relativity in a regime which 
is well beyond our current technological capabilities for detection. Henceforth, the emphasis of these 
studies is not as much in predicting new phenomena that could be tested in a short period of time as 
it is in gaining insight on how gravitational lensing in particular, and gravity in general, behave in the 
strong field regimes. 

In ordinary lensing situations, the lens is placed between the source and the observer. But if the 
lens is a compact object with a photon sphere and the observer is placed between the source and 
the lens, it leads to the formation of images with deflection angles closer to odd multiples of vr. The 
observer sees the images in front and the source behind. This situation is called retro-lensing and was 
studied for the first time by Holtz & Wheeler jj^, who analyzed only the two stronger images for a 
black hole placed in the galactic bulge with the sun as source. They also proposed retro-lensing as a 
new mechanism for searching black holes. A similar idea was suggested in the context of defocusing 
gravitational lensing by Capozziello et al JH]- De Paolis et al Jlj recently analyzed the retro-lensing 
scenario for a bright star close to the massive black hole at the galactic center. 

In this paper we give a complete treatment of retro-lensing by a static, spherically symmetric 
lens, using the strong field limit. In Section 2 we obtain the lens equation for retro-lensing and the 
defiection angle in the strong field limit. In Section 3, the position and magnification of the images are 
calculated. In Section 4, the results are compared with those for ordinary lensing and astrophysical 
examples are given in Section 5. Final concluding remarks are given in Section 6. Throughout the 
paper we use units such that G = c = 1. 

2 Lens equation and deflection angle 

Consider an observer (o) placed between a point source of light (s) and a strong field object with a 
photon sphere (e.g. a black hole), which we will call the lens (1). The line joining the observer and 
the lens define the optical axis. The background space-time is considered asymptotically flat, with 
the observer and the source immersed in the flat region. We call f3 the angular position of the source 
and 6 the angular position of the images (i) as seen by the observer. The lens situation is shown in 
Fig. ^ We can take /3 > without loosing generality. The lens equation for retro-lensing is slightly 
different from ordinary lensing (i.e. when the lens lies between the source and the observer), due to 
the different relative positions of the lens, source and observer: 



where dos = Dos/2M and dig = Dis/2M are, respectively, the observer-source and the lens-source 
distances in units of the Schwarzschild radius, and a is the deflection angle. In lensing situations where 
the objects are highly aligned, the angles /? and 6 are small and a is closer to an odd multiple of tt. For 
the images at the opposite side of the source (see Fig. P), it can be written as a = (2n — 1)-k + Aa„, 
with n G N and < Aa„ <C 1. Then, the lens equation takes the form 




(1) 




(2) 
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Figure 1: Lens diagram. The observer (o), the lens (1), the source (s) and the projection on the lens 
plane of the first relativistic image (i) are shown. The angular position of the source is /3, of the 
image is 9, and doi, dos, dig are, respectively, the observer-lens, the observer-source and the lens-source 
distances. Two infinite sets of relativistic images (with both directions of rotation around the lens) 
are obtained with deflection angles close to odd multiples of vr. 

As in the case of ordinary strong field lensing, two infinite sets of relativistic images are formed. To 
obtain the other set of images (at the same side of the source) we should take a = — (2n — l)7r — Aa^, 
so Aa^ must be replaced by — Aa„ in Eq. In the case of perfect alignment, an infinite series of 
concentric Einstein rings are obtained. 

To calculate the defiection angle, consider an static spherically symmetric metric 

ds^ = -f{x)dt^ + g{x)dx^ + h{x)dn^, (3) 

where x = r/2M is the radial coordinate in units of the Schwarzschild radius. We require this metric 
to have a photon sphere, which radius Xps is given by the greater positive solution of the equation: 

f}M = iM. (A) 

h{x) fix) ■ ^ > 
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where the prime means derivative with respect to x. 



The deflection angle corresponding to the images situated at the opposite side of the source as a 
function of the closest distance of approach xq is given by jT^ 

a{xo) = /(xo) - IT, (5) 

with 

The impact parameter b (in units of the Schwarzschild radius) is ^S] 

Prom the lens geometry we also have that 

b{xo) = doi sin ^ doiO, (8) 
with doi = Dqi/2M the observer-lens distance in units of the Schwarzschild radius. 

Bozza has shown that for a spherically symmetric lens with a photon sphere, light rays with 
impact parameter h close to the photon sphere have a deflection angle that can be approximated by 

a = -ai\n{b-bps) + a2 + 0{b-hps), (9) 

where bps = b{xps), and ai^2 are constants which depend only on the type of lens. Provided the metric 
of the compact object that acts as a lens, the coefficients ai^2 can be obtained in a laborious but 
straightforward way using the formalism developed in [2]. The method consists essentially in sepa- 
rating the integral of Eq. © in a divergent part which gives the logarithmic term in Eq. Q, and a 
regular part which gives the constant a2- This approximation, named the strong field limit, gives very 
accurate results. Por more details see 

Another way of calculating the constants ai^2 is following a similar approach to that used by Eiroa, 
Romero & Torres [S]: 

a,= hm -[H^o)-bM]r{xo) ^ ^^^^ 
xo^xps b'{xo) 

and 

as = In I lim \e'^^^^-^[b{xo) - b{xps)r]\ , (H) 

with the prime now meaning the derivative with respect to xq. These limits can be evaluated numer- 
ically with standard software. 

The values of the coefficients oi and a2 for Schwarzschild and Reissner-Nordstrom lens geometries 
are given in the Appendix. 
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3 Position and magnification of the images 

Inverting Eq. @ to obtain the impact parameter 

' 02 — a 



b = bps + exp 



ai 



and using that b = doi6, we have that the position of any image is 

bps I / 02 - a 

6* = + — exp 

doi doi V «i 



(12) 



(13) 



Taylor expanding to first order around a = (2n — 1)tt, we can approximate the angular position of the 
n-th image by 



lO A 



with 



and 



- doi + doi P 



02 — (2n — l)7r 



= — exp 

Ol 



Ol 

02 — (2n — l)7r 

Ol 



Then from Eq. (|14jl we obtain 



i^Oln = 7 doU 

sn 



and replacing it in the lens equation 



dosdoi I Q, 

aisKn 



The last equation can be written in the form 



1 + 



dlsS,n 



and as the second term in the parentheses is much greater than one, we have that 

dlst,n 

We finally obtain that the angular positions of the images are 



Jn — 



dn 



(14) 
(15) 
(16) 

(17) 
(18) 

(19) 

(20) 

(21) 



The second term in Eq. (|21jl is a small correction on 6^, so all images lie very close to 6^. With a 
similar treatment, the other set of relativistic images have positions given by 



(i'OsO'ol 

When /3 = an infinite sequence of Einstein rings is formed, with angular radius 

aE nO disin 
^os^ol 



(22) 



(23) 
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Gravitational lensing conserves surface brightness, so the ratio of the sohd angles subtended by the 
image and the source gives the amplification of the n-th image: 



Using Eq. (|21|) we have that 



sin (3 dj3 


-1 


/3 df3 


-1 


sin On dOn 




On dOn 





1 



,0 ^ MlL(^p _ 00) 



dos doi 



dos doi 



which can be approximated to first order in dis/dosdoi <C 1 by 

_ 1 dis 



(24) 



(25) 



(26) 



(3 dosdoi 

The same result is obtained for the other set of images. The first relativistic image is the brightest 
one, and the magnifications decrease exponentially with n. 

oo 

The total magnification, considering both sets of images, is /i = 2 ^ which can be written as 

n=l 



1^ 



2 distts 



(3 dosdl,ai ' 



where 



as = h 



exp 



a2— TT 
ai 



ps 



+ 



exp 



0.1 



(27) 



(28) 



Extended sources 



If the source is extended, we have to integrate over its luminosity profile to obtain the magnification: 



(29) 



where I is the surface intensity distribution of the source and is the magnification corresponding 
to each point of the source. If the source is uniform, and we use polar coordinates (ii, ip) in the source 
plane, with i? = in the optical axis, and we consider the source as a disk D{Rc, Rs) of radius Rs 
centered in Rc, the last equation takes the form 



and, using that (5 = R/Dqs is the angular position of each point of the source, we have that 

//d(/3,,a) ^^pPdpd^ 



where D{Pc,Ps) is the disk with angular radius (3s centered in f3c- Then 

fJ-n 



I dls r,Oc 



■Kj3l dosdoi 



(30) 



(31) 



(32) 
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and 



2/ distts 



with / = Jjj^ii^ p \ dpdip. This integral can be calculated in terms of elliptic functions 1^: 



/ = 2Sign[ps - f3c 



2'(/3s-/3c)2 



where 



r4>o _i 
F(0o,A)= / (l-Asin2(/>) 2 
Jo 



E{(t>o,X) 



'I - Asin^ 



(33) 

(34) 

(35) 
(36) 



are respectively, elliptic integrals of the the first and second kind with the arguments (pQ and A. 



4 Comparison with ordinary lensing 

In ordinary lensing the angular positions of the two sets of relativistic images, for a point source, are 
given by 

+ ^{P-e% (37) 



^ _ ^0 _|_ dps^n 

disdoi 



and 



with 



and 



n = -0° + ^(/3 + O, 

1 / a2 — 2mT 



90^^ + exp 

" doi doi 



ai 



1 f a2- 2n7r\ 
U = — exp ( I 



(38) 



(39) 



(40) 



ai \ oi 

We see that for both cases (ordinary and retro-lensing) , the images are formed close to 9^, which has 
values of the order of the angular radius of the photon sphere of the lens. From Eqs. (|15() and ()39|) we 
observe that the images for retro-lensing lies a bit farther of the photon sphere than in the ordinary 
case. 



The magnifications in ordinary lensing, for a point source, of the relativistic images are (both sets): 



and the total magnification, summing the magnifications of both sets of relativistic images, is 

2 dosO-3 

(5 disdlfli 

where 

(43, 



1-exp -f 1-exp -f 
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Comparing Eqs. (|28|) and ()4c{() we see that the first term of 03 in retro- lensing is greater by a factor 
g7r/ai g^j^i^ ^j^g second by a factor e^'^/'^^ with respect to the ordinary case. Then, taking by example 
the lensing by a Schwarzschild black hole, and considering a situation where the factor involving the 
distances and the source angular position is the same, the magnification of the relativistic images in 
retro-lensing is greater by a factor about 24 compared with ordinary lensing. 

In addition, with ordinary lensing, besides the two set of relativistic images, two weak field images 
with small deflection angles are produced by the lens. These weak field images do not appear in retro- 
lensing, and it is the main difference in both cases. The weak field images are orders of magnitude 
stronger than the relativistic ones [^IS] and lie very close to them. This makes detecting the relativistic 
images in ordinary lensing a very difficult task. In the case of retro-lensing, instead, the faint relativistic 
images are stronger and relatively easier to identify due to the absence of the brighter weak field images. 

5 Examples 

In this Section we apply the formalism to the Schwarzschild and Reissner-Nordstrom lens geometries. 
Two cases are analyzed: a super-massive black hole situated at the galactic center and a stellar mass 
black hole in the galactic halo. In both cases a nearby star is taken as source. 

There is strong evidence of the existence of a super-massive black hole at the center of our galaxy 
[ini. We take a M ~ 2.8 x 10^ Mq black hole as an fiducial lens at a distance Dqi = 8.5 kpc. A 
star with radius R = Rq situated at Dqs = 50 pc is assumed as the source. In Table ^ the first 
(outer) Einstein ring angular position {9f) and the limiting angular value of the Einstein rings {0^) 
are given for the Schwarzschild and Reissner-Nordstrom black hole lenses. The images lie very close to 
the Einstein angular radii. For extended sources, when there is complete alignment, instead of a pair 
of images for each n (one on each side), an annular shaped image is formed, with magnification 2^^. 
The magnification of the first relativistic image and the total magnification, for perfect alignment, are 
also shown in Table ^ . 



Table 1: Black hole in the galactic center as a gravitational retro-lens for a star situated at 50 pc from 
the Earth. For several values of the charge q (in units of the Schwarzschild radius), the angular radii 
of the photon sphere {9ps), the first Einstein ring (^f"), and the limiting value of the Einstein rings 
(0^), are given. The magnifications, for perfect alignment and taking the source as extended with 
radius Rq, of the first image is 2/ii, whereas the total magnification is /i. 



kl 





0.1 


0.2 


0.3 


0.4 


0.5 


0ps(^arcsec) 


9.755 


9.667 


9.395 


8.899 


8.080 


6.503 


(//arcsec) 


17.39 


17.28 


16.95 


16.37 


15.53 


15.04 


^^(^arcsec) 


16.90 


16.78 


16.43 


15.80 


14.78 


13.01 


2iUi(xlO-i°) 


3.032 


3.031 


3.049 


3.164 


3.673 


7.705 


/x(xlO-iO) 


3.037 


3.037 


3.055 


3.172 


3.687 


7.784 



Low mass black holes are common in our galaxy Let us now consider a stellar mass black 

hole with mass M = 7Mq in the galactic halo with Dqi = 4 kpc as retro-lens and a star situated at 
Dos = 50 pc with radius R = Rq as source. The results are given in Tabled 
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Table 2: Same as Tabled for a black hole in the galactic halo as a gravitational retro-lens with a star 
situated at 50 pc from the Earth as source. 



kl 





0.1 


0.2 


0.3 


0.4 


0.5 


^ps(x 10~^/xarcsec) 


5.182 


5.136 


4.991 


4.728 


4.292 


3.455 


6'f (xlO^^/iarcsec) 


9.236 


9.178 


9.002 


8.696 


8.249 


7.990 


^^(x 10~^^arcsec) 


8.976 


8.916 


8.729 


8.393 


7.853 


6.910 


2^i(xl0-2i) 


4.053 


4.052 


4.076 


4.230 


4.911 


10.30 


/i(xl0-2i) 


4.060 


4.060 


4.084 


4.240 


4.928 


10.41 



In both examples, the dependence of the angular positions and magnifications of the images with 
the charge is weak, obtaining important changes only for extreme values of charge. 

As the images lie close to the photon sphere, the angular positions of them take very small values, 
what increase with larger lens mass and smaller observer-lens distance. For the supermassive black 
hole at the galactic center, they are of the order of 10 //arcsec, and for the stellar mass black hole in the 
galactic halo they are about 10~^ yuarcsec. Something similar happens with the magnifications, which 
also grow with greater lens mass and smaller observer-lens distance. For the supermassive black hole 
at the galactic center the magnifications are of order 10^^^ and for the low mass one in the galactic 
halo of order 10^^^ . 

As long as dos ^ doi, the angular positions and amplifications of the images are not very sensitive 
to the value of dos- So if we take the Sun or a star at 100 pc from the Earth as source, the order 
of magnitude of the results shown in Tables H and [21 does not change. As the images appear highly 
demagnified, to improve the possibility of observing them, strong nearby sources with nearly perfect 
alignment are needed. 

In all cases, the first image is the strongest one, carrying about 99% of the total luminosity. 

6 Concluding Remarks 

We have presented the strong field limit formalism for gravitational retrolensing. We have found 
the positions and magnifications of the relativistic images within the strong field limit, and studied 
some cases which might present particular interest and be a benchmark point for comparison with the 
normal lensing situation. As in the case of ordinary lensing, the relativistic images are faint. However, 
if we compare the magnification of retro-lensed relativistic images with those of the ordinary strong 
field lensing, we find that the former are significanlty stronger. Additionally, as no brighter weak 
field images appear in retro-lensing, it would be easier to detect them. The angular separation of the 
relativistic images is of course beyond the angular resolution of current or foreseeable technologies, 
and not only angular resolution but also sensitivity improvements are needed to detect them (see the 
discussion in Eiroa, Romero & Torres jSj). The image separations and magnifications present values 
that make the observational identification of these relativistic images a challenge for the future: we are 
not aware of any current research, in any wavelength, which would allow to resolve the retro-lensing 
images in the next years. In the case of charged black holes, it is only for extreme values of charge 
that important changes are obtained. 
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Appendix: Values of the coefficients ai and a2 

Schwarzschild black hole 
The metric functions are 

f{x)=g-\x) = l--, h{x)=x\ (44) 

X 

then 

3 

Xps = ^, (45) 



I ) bps — 

xoj '^2 



b{xo) = xo ( 1 - — ) , bps = — (46) 



The coefficients can be calculated exactly |1J [S] 

ai = 1, (47) 



02 = ln[324(7\/3 - 12)] - tt. 



(48) 



Reissner-Nordstrom black hole 



1 



In this case 

f{x)=g~^{x) = l-- + ^, h{x) = x^, 

X 

where q = Q/2M is the electric charge in units of the Schwarzschild radius. Then 



X. 



ps 



32 9 



1/2- 



1 

6(xo) = xo ( 1 \- —r, 



b. 



(3 + ^9 - 32g2)2 



'ps 



4^/2^3 - 8g2 + y/g - 32g2 
The coefficient oi can be calculated exactly 



ai 



(3 - Xps)xjg - dq^Xps + 8g4 



(49) 

(50) 

(51) 
(52) 



(53) 
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and the coefficient 02 can be approximated [Hj by 

02 = ci + C2 + ai In bps - vr (54) 

with 

ci = 21n[6(2-^/3)] + ^{^/3-4 + ln[6(2-^/3)]}<7' + 0(/), (55) 



C2 = ai In 



2{xps - q 



{xps - 2q^)^{xjs - Xps + g2) 



(56) 
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